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Abstract 

We study the cosmology of the Supersymmetric Standard Model augmented by a 
gauge singlet to solve the //-problem and describe the evolution of the domain walls 
which are created during electroweak symmetry breaking due to the discrete Z% 
symmetry in this model. The usual assumption, that non-renormalizable terms in- 
duced by gravity (which explicitly break this symmetry) may cause the walls to 
collapse on a cosmologically safe timescale, is reconsidered. Such terms are con- 
strained by considerations of primordial nucleosynthesis, and also by the fact that 
by not respecting the Z% symmetry they induce divergences which destabilise the 
hierarchy and reintroduce the //-problem. We find that, even when the Kahler po- 
tential is 'non-minimal' (i.e. when the hidden sector couples directly to the visible) 
the model is either ruled out cosmologically or suffers from a naturalness problem. 
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1 Introduction 



The purpose of introducing (softly broken) supersymmetry into the Standard Model is 
to bring under control the quadratic divergences associated with a fundamental Higgs 
boson and make it 'natural' for its mass to be at the electroweak scale M. Yet the 
minimal supersymmetric Standard Model (MSSM) has its own naturalness problem. Its 
Lagrangian contains a term // H\H 2 mixing the two Higgs doublets which are now required 
to give masses separately to the up- and down- type quarks. For successful phenomenology 
\x should also be of order the electroweak scale but this must now be set by hand — the '/i- 
problem' pi, |3|]. To address this problem, the next-to-minimal supersymmetric Standard 
Model (NMSSM) Q contains an additional singlet Higgs superfield N. By invoking a Z 3 
symmetry under which every chiral superfield $ transforms as $ — > e 27n / 3( 3>, the allowed 
terms in the superpotential are now XNHiH 2 — |iV 3 (in addition to the usual Yukawa 
terms generating fermion masses) while the Higgs part of the soft supersymmetry breaking 
potential is extended by the inclusion of two additional trilinear soft terms A\ and to 

V^ s = -\A x (NH 1 H a + h.c)-£A k (N*+h.c.) 

+m 2 Hi \H 1 \ 2 + mj l2 \H 2 \ 2 + m 2 N \N\ 2 , (1) 

where H\H 2 = H^H® — H~H + . The //-term can now be simply set to zero by invoking 
the Z 3 symmetry. An effective /x-term of the form X(N) will still be generated during 
SU(2)xU(l) breaking but it is straightforward to arrange that (N) is of order a soft 
supersymmetry breaking mass. Apart from solving the '/i-problem' the NMSSM also has 
interesting implications for supersymmetric phenomenology |5j and dark matter [[J. 

Unfortunately, the NMSSM runs into a cosmological difficulty The Z3 of the model 
is broken during the phase transition associated with electroweak symmetry breaking in 
the early universe. Due to the existence of causal horizons in an evolving universe, such 
spontaneously broken discrete symmetries lead to the formation of domains of different 
degenerate vacua separated by domain walls |7|, p. These have a surface energy density 
o ~ v 3 where v is a typical vacuum expectation value (vev) of the fields, here the elec- 
troweak scale of C(10 2 ) GeV. Such walls would come to dominate the energy density of 
the universe and create unacceptably large anisotropies in the cosmic microwave back- 
ground radiation unless their energy scale is less than a few MeV Therefore cosmology 
requires the Z 3 walls to disappear well before the present era. Following the original sug- 
gestion by Zel'dovich et alM, this may be achieved by breaking the degeneracy of the 
vacua, eventually leading to the dominance of the true vacuum. This happens when the 
pressure, i.e. the difference in energy density between the distinct vacua, begins to exceed 
the tension a/R, where a is the surface energy density of the walls and R the scale of their 
curvature. When R becomes large enough for the pressure term to dominate, the domain 
corresponding to the true vacuum begins to expand into the domains of false vacuum and 
eventually fills all of space. It was recently argued |L0| that gravitational interactions at 



the Planck scale Mpi would explicitly violate any discrete symmetry, causing just such 
a non-degeneracy in the minima of Oiy' /Mpi) where v is a generic vev (of 0(Mw) in 
our example). In fact, this suggestion had been applied already to the NMSSM in the 



context of string theories [11]. Thus there would appear to be a natural solution to the 



cosmological domain wall problem for the NMSSM. 
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In this paper we study whether this solution is indeed viable. In the following section 
we derive the structure of the walls, and show that the surface energy is approximately 
as expected on dimensional grounds. We go on to describe the evolution of the walls 
under the influence of the tension, the pressure due to the small explicit Z 3 breaking and 
the friction due to particle reflections. In particular we demonstrate that wall domination 
of the energy density of the universe is avoided if the gravitationally induced terms are 
of order six or less. This is not however the tightest constraint on the domain walls; by 
applying constraints based on primordial nucleosynthesis we show that the magnitude 
of Z 3 breaking must be £ 10~VM^/Mpi, in order to make the walls disappear before 
the nucleosynthesis era begining at T ~ 1 MeV. Thus only operators of dimension five 
(suppressed by at most one power of the Planck mass) are permitted. This reduces to 
three the number of possible Z 3 breaking terms which are allowed in the superpotential 
or Kahler potential and which induce dimension-5 operators in the effective potential. By 
inspection we find that the existence of one or more of these operators implies that there is 
no symmetry (discrete, global, gauged, i?-symmetry or gauged .R-symmetry) under which 
the low-energy singlet can be charged. This implies that there cannot be any explanation 
for the absence of three allowed low energy operators which include the /z-term itself 
as well as quadratic and linear terms in N. Thus our first conclusion is that not only 
does the NMSSM not solve the /i-problem, it actually makes things worse by introducing 
additional operators and by disallowing any symmetry which would forbid them. 

We then go on to consider the fact that the singlet which appears in the NMSSM may 



introduce destabilising divergences [12]. Essentially the problem is that by introducing 
non-renormalizable terms together with soft supersymmetry breaking, we may produce 
corrections to the potential which are quadratically divergent and thus proportional to 
powers of the cut-off A in the effective supergravity theory. Since the natural scale for this 
cut-off is Mpi, these can in principle destabilise the hierarchy, forcing the singlet vev and 
hence the scale of electroweak breaking to become very large (at least of order y/M^Mpi) . 
By examining the possible Z 3 breaking terms, we demonstrate that the removal of domain 
walls by this mechanism indeed destabilises the hierarchy. We conclude that the two 
constraints, viz. stability of the hierarchy and domain walls, cannot be simultaneously 
satisfied by any gravitationally suppressed operators which one can add to the Lagrangian. 

We consider alternative ways for dealing with the domain walls. One possible solution 
is to reintroduce the \i term in the superpotential in such a way as to avoid the introduction 
of the dangerous non-renormalisable operators. If we drop the assumption of minimality 
in the Kahler potential by allowing certain couplings of the hidden sector fields to the 
visible sector (as in Ref.||) we can retain Z% symmetry in the full theory, but break it 
spontaneously when supersymmetry is broken. In this way the hierarchy is not destabilised 
by tadpole diagrams. However the naturalness problem cannot be solved even for these 
more general models. 

Finally we consider how gauge singlets may be accommodated in supersymmetry, 
without invoking these problems. There appear to be only a few possibilities, none of 
which yields a phenomenology bearing any resemblance to the NMSSM. 
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2 Domain Walls in the NMSSM 



When a discrete symmetry is spontaneously broken as the universe expands and cools, 
'domains' of the different degenerate vacua form separated by narrow regions of higher po- 
tential called 'domain walls' ||. The structure of these walls may be determined by finding 
time-independent solutions to the classical field equations after imposing the boundary 
conditions that at the endpoints the fields should be in distinct vacuum configurations. 
This has been done using numerical methods for the NMSSM potential |13] and we re- 



iterate the essential features of the Z 3 walls. As might be expected from dimensional 
arguments and by analogy with the analytically soluble case of a single real scalar field 
in a Z 2 symmetric potential |§, the thickness and energy density of the walls are of order 
v~ x and z/ 3 respectively, where v is a typical vacuum expectation value. For naturalness 
reasons one would tend to assume that all three vacuum expectation values are of the 
same order; however, it is also possible that the singlet vev, x, is much larger than the 



usual v = y isf + i/f — 174 GeV. This is in fact quite likely in the light of recent analyses 
where unification of soft terms and gauge couplings is imposed at the GUT scale; the only 
viable scenarios are then found to have xjv >, 10 with especially large values when the 



Higgs sector Yukawa couplings are very small [fl4[ . In such cases, we would expect the 
wall to have a much higher surface energy a; indeed we find that this is well approximated 
by 

a~5xl0 7 GeV 3 (-^) ( — V , (2) 



v 0.1 J \5u, 

when x is at least a few times larger than v. This formula is accurate to about a factor 
of 2 in practice and is very good for large x, relative to both the trilinear soft terms and 
to v. 

Similarly the thickness of walls is given by 

^2xlO-GeV-> (A)" (»)-' (3 ) 

which again is most accurate when x 3> v and x 3> Ak, A\. We show an example of a 
wall with large x in Figure 1. In comparison to the cases shown in Ref.||13||, we see that 
the wall is thinner and the surface energy higher, as expected. 

We note that if both Ak and A\ are zero, then the Z% symmetry of the scalar potential 
becomes a U(l) symmetry, and so the wall energy falls to zero and its width becomes 
infinite; in this limit however we have an axion problem. We find that if A^ or A\ are 
greater than a few GeV then the wall energy is insensiteve to their exact values. 

Immediately after the electroweak phase transition the universe is filled with equal 
volumes of the three degenerate phases. These are correlated on a length scale which 
depend on the nature of the phase transition, varying from £ ~ T^ 1 for a second-order 
transition to £ ~ H" 1 for a strongly first-order transition ||, [TjJ. Since the probability 
for each vacuum (0.333) is just above the percolation threshold (which for continuum 
percolation theories is found to be 0.295 [16|]), the universe is then filled with highly 
convoluted, infinite regions separated by stable domain walls of typical curvature scale £, 
which rapidly grows to the size of the horizon. 

Let us now turn to the dynamics of cosmological networks of such walls. As discussed 
in Ref. M, the most important forces acting on the walls are surface tension, friction and 



4 



pressure. The equation of motion for a quasi-spherical piece of wall moving with velocity 



(4) 



R (with 7 = 1/y 1 — R?) and having local radius of curvature R, is 

d 2 R 2 (npv) e 
dt 2 R'-f 2 o"7 3 o"7 3 



The first term on the RHS reflects the fact that it is energetically favourable for the 
wall network to reduce its surface area through surface tension, and hence small domains 
will collapse, irregularities in the surfaces will straighten out, and the correlation length 
will increase. This term expresses just the conservation of energy in the absence of pressure 
and friction. 

The second term on the RHS corresponds to friction arising due to the interactions of 
the wall network with the thermal plasma. As particles reflect off the walls, they exert a 
force given by the thermally averaged momentum transfer (npv), where n is the particle 
density, v the particle velocity relative to the wall, and p the momentum perpendicular 
to the wall. (Actually the friction is oc v only when »<c.) Friction is clearly important 
at times very close to the electroweak phase transition if the top quark and gauge boson 
fields are still in equilibrium in the plasma. At later times, when the number density of 
these particles is exponentially suppressed, the main source of friction is the interaction 
of the walls with lighter fermions in the plasma. The constant difference in phase in 
the mass terms on either side of the wall (i.e. ir/3 or 2n/3) does not by itself cause 
any reflection but rather just a phase shift in the fermion masses (as can be checked by 
equating transmission and reflection coefficients at the wall). In order to estimate the 
reflection coefficient, it is useful to describe the space dependent mass by the inverted 
bell-shaped function 

m 2 (x±) = m 2 — a 2 ^\ — , (5) 

cosh ax± 

where x± is the perpendicular distance from the wall, and m is the mass given to the 
reflecting particle by the Higgs fields which comprise the domain wall of width a -1 . The 
task of finding the reflection coefficient (using the Klein-Gordon equation) then reduces to 
a known problem, the modified Poschl- Teller potential, which can be solved analytically 
(see for example Ref.[17j). We take the depth of the well to be m 2 and the width a -1 ~ 



My/. The depth parameter is A = (1 + m 2 /M w ) and the reflection coefficient is then 
found to be 

2 4 4 
\J\\ l — ~ (Q) 

1 1 vr 2 m 4 + M w sinh 2 irpa ~ p 2 M w ' v ; 

where we have taken Mw P 3> m as is appropriate once the gauge bosons and top 
quark have fallen out of equilibrium. (There is a region at low energy \p\ < m 2 /Mw in 
which the particles experience total reflection ||. However this contribution is insignif- 
icant here, being suppressed by many powers of m 2 /Mw-) Clearly particles which are 
heavy, especially the bottom quark, will be more important here. We can estimate the 
friction by considering a particle of mass m, when the wall velocity through the plasma, 
u, is small. Then 



(npv) = g 



d 3 p m 4 p\ 
{2itYp 2 M w ~E 



1 1 



exp^E + 7wp_|_)/T) + 1 exp^E — , yup±)/T) + 1 



(7) 



where T is the temperature of the plasma, and g is the number of degrees of freedom of 
the reflecting particles. Expanding this function in u keeping the leading term only and 
performing the angular integral, we find 



where 



F(x m ) — I dx x m (x x m ) 2 , (9) 



and we have defined x = E^jT and x m = mj/T. This integral is very well approximated 
by 

F(x m ) =a4(0.6e-^) 3 . (10) 
Summing over all the particle species in the plasma, we find that 



(npv) = f(T) u ^, (11) 



where f(T) < 5 x 10~ 4 at all temperatures. We show f(T) in Figure 2 where, apart from 
omitting the contribution of the up and down quarks, we have neglected the possible effect 
of the quark-hadron phase transition. In the era when pressure is negligible (i.e. when 
the typical curvature scale is small), we can calculate the terminal wall velocity, Mterm, 
and establish a posteriori that our approximation of small u to obtain eq. (Q) is indeed 
correct, i.e. friction is important. Substituting the friction into eq.(f|), we find 

16tt 2 / a \ 

Mterm = • (12) 

For typical values of the radius, R ~ Uterm^, we see that friction is important only at 
temperatures above a few hundred MeV. We therefore conclude that shortly after the 
quark-hadron phase transition the walls move with velocities comparable to the speed of 
light and so we may safely neglect friction in what follows. 

The last term on the RHS in eq.(|]) is the pressure corresponding to the difference e 
in the energy density between the different vacua. As remarked earlier, this will become 
dominant when it exceeds the surface tension, i.e. when 

e>^- (13) 

We show this happening in Figure 3, where we have performed a simple thin wall sim- 
ulation of a network of domain walls using techniques similar to those used in ref . fl8fl . 



and which we have discussed in more detail elsewhere |13| . In the absence of friction it is 



convenient to rescale the parameters with some typical length scale, R , which we choose 
to be 1 cm, corresponding approximately to the curvature scale when pressure becomes 
dominant if e ~ M^/Mpi. Thus defining p = R/Ro, and r = t/Ro, eq.(§) becomes 

d 2 p = 2_ _ £^o , u x 

dr 2 fry 2 7% 
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Thus there are only two independent parameters in our simulation, given by the pressure 
in each of the two false vacua, eRo/a. Initially, the walls expand under their own tension, 
and the structure develops in the manner discussed in Refs. |L8|, [19| . Eventually pressure 
dominates as expected and the entire volume is cleared of walls. This contrasts with the 



no-pressure case, where one or two horizon-sized walls always remain [13]. The behaviour 



for different values of the pressure or surface tension is identical if the time and length, 
respectively, are scaled appropriately. 

One might consider the possibility that since frictionless, pressureless walls expand 
until there is roughly one wall per horizon scale |§, domain walls may be accomodated by 
simply assuming that our local region of space-time just happens to be empty of them, 
i.e. that there is a wall lurking just outside our present horizon. There are at least two 
objections to this. Firstly the walls eventually come to dominate the energy density of the 
universe, causing unacceptable 'power- law' inflation | 20fl , unless their separation is many 



times greater than the present horizon scale, which is clearly impossible by causality. 
Secondly, even such a wall outside the horizon will have a curvature scale comparable to 
the present horizon scale and thus induce unacceptably large anisotropy in the cosmic 
microwave background ET . 



3 When Walls Collide 

What value of the pressure (i.e. explicit Z 3 breaking) is required to safely remove the 
walls? The crudest estimate we can make is simply to insist that the walls are removed 
before they dominate over the radiation energy density in the universe in order to avoid 
wall driven inflation. Since the walls move at close to the speed of light below the quark 
hadron phase transition, their curvature scale will be roughly the horizon size, R ~ 
t~ M Pl /gl /2 T 2 . Since the energy density of the walls is 

a 

Pwalls ~ , (15) 

and the radiation energy density is p ra d ~ gl^ 2 T 4 , we see that walls dominate the evolution 
below a temperature 

/ \ 1/2 

T* ~ • (16) 

To prevent this we require the pressure to have become dominant before this epoch, i.e. 

a a 2 , , 

£> 5T~^ (17) 

A pressure of this magnitude would be produced by dimension-6 operators in the potential. 
However, one should consider further constraints coming from primordial nucleosynthesis, 
and we find that only operators of dimension-5 or less are sufficient to satisfy these. In fact 
for weak scale walls the time associated with the temperature is i* ~ M^>J gV 2 ~ 
10 8 sec, i.e. long after nucleosynthesis. The entropy produced when the walls collide 
(which is by now a major proportion of the total entropy in the universe) is dumped into 
all the decay products of neutral Higgs particles, i.e. Standard Model quarks and leptons. 
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In order to check whether this violates phenomenological bounds, we compute the relative 
energy density released in such collisions, viz. 

Pwalls ° 



7 x 10" 11 GeV 



a 



Mir, 



sec. 



where we have taken the number of relativistic degrees of freedom in the plasma to be 
g* = 43/4. Detailed consideration of the effects of high energy particles on primordial 
nucleosynthesis and on the 2.73 K Planckian spectrum of the microwave background ra- 
diation impose severe upper limits on this parameter For the typical values of o in 
Figure 1, we find that the walls are required to disappear before the onset of nucleosyn- 
thesis at about 0.1 sec, as otherwise the hadrons in the showers triggered by the decay 
products would alter the neutron-to-proton ratio, resulting in a 4 He mass fraction in ex- 
cess of the conservative observational upper bound of 25% |23[]. This means that in order 
not to disrupt primordial nucleosynthesis, we require explicit Z% breaking of magnitude 

e ~ AVM^/M p1 , (19) 

with 

A'>1(T 7 . (20) 



4 The Return of the fi Problem 

Having established that one needs dimension-5, Z 3 breaking operators to appear in the 
effective potential, we can consider ways in which this can be achieved by adding terms 
to the Kahler potential or superpotential. We first assume that these are 'minimal' in the 
sense that they do not contain couplings between the hidden and visible sectors (which 
couple only through gravity). Later we shall consider the most general non-minimal 
case. In all cases we find that there is a naturalness problem associated with the explicit 
breaking of the Z 3 symmetry. 

Let us write down the contributions to the supergravity Lagrangian which explicitly 
break Z 3 , and which are invariant under the NMSSM gauge group. These are 

A Mpj' A M P1 ' A ^WpT' (21) 



in the superpotential, and 

+- \'')(H;H' ! ' ) 



ai ' H )/' (22) 



in the Kahler potential. As in Ref. ||12||, we can absorb the last two contributions into the 
superpotential to O(Mp^) by making the redefinitions 

V M Pl J 

N - iV-^M, (23) 
M P1 ' 1 ; 
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and so we shall consider only the first three contributions in what follows. Inspecting 
these, we observe that N must be a singlet under any additional symmetry in order 
for any one of these terms to exist in addition to the terms N 3 and NH1H2 in the low 
energy superpotential. In other words, each of them implies that the following 'unnatural' 
contribution to the superpotential is invariant 

^'unnatural' = »" N + ^ + A*^- ( 24 ) 

Thus not only have we reintroduced the /i-problem, we now have two additional natural- 
ness problems. Whereas the standard /^-problem may well be solved at a future date (for 
example by the mechanism of Ref. 0), we shall see that the naturalness problem which 
has reappeared here can have no solution based on an underlying symmetry. 



5 The Return of the Hierarchy Problem 

As if the difficulties above were not bad enough, there is the possibility of quadratic tadpole 
divergences which can lead to a destabilisation of the hierarchy [|12|]. This exacerbates 
our problems, since such divergences arise at each order in perturbation theory, forcing 
us to re-fine-tune. These are a potential problem in any supergravity model with gauge 
singlets since the dangerous diagrams are not excluded by gauge invariance. 

These diagrams arise when supersymmetry is spontaneously broken, because super- 
Weyl-Kahler invariance necessitates that the vev of the Kahler potential become non- 



trivial. In fact 12 



(e 2 */ 3 ) « e 2K / 3 \(l + 9 2 M 2 + 9 2 M 2 + 6 2 9 2 M*) 

(0) « 0|(l + fl 2 M s 2 ), (25) 

where is the chiral compensator, M s is the scale of supersymmetry breaking in the 
hidden sector, and the RHS refers to only the scalar components. The leading tadpole 
divergences are quadratic and appear at two-loop order for the first two operators in 
eq . (|2l|) . In our case, the diagrams responsible are shown in Figures 4a and 4b, and they 
lead to the terms 

^^y^N + <P* N )M Pl m 2 3/2 + ^^(F N + F* N )M m m 3/2 (26) 

and 

A'A , , o A'A 



( , () -> )2 + <P* N )Mpimi /2 + -^y(F N + F* N )M Pl m 3/2 (27) 
respectively, where we have taken the cut-off to be A ~ Mp\ and introduced the gravitino 



mass m 3 / 2 ~ yM s 2 /Mpi. Here A and k are the Higgs sector Yukawa couplings defined 
earlier. The third term in eq.(pTf) gives rise to a divergence at three- loop order as shown 
in Figure 4c and the calculation is a little more tricky. Using the perturbation theory 
rules of Ref . [ T2j , quadratic divergences are indeed found to arise of the form 



^^ 3 (<pN + <P* N )M Pl ml /2 {21 
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and 

j0^(F N + F* N )M PimV 2, (29) 

where we have replaced a quadratically divergent three loop integral with a cut-off, Mp V 
All of these terms naturally drive the vev of the singlet (and hence of Hi, H 2 ) to the 
hidden sector scale, (x) ~ Jmz/2M-p\ ~ 10 11 GeV. If we wish to avoid the reappearance 



of the hierarchy problem, these terms should be smaller than ~ (0jy + 0Af) m 3/2 or ~ 
(Fjv + F^)m^ 2 . Even for the three loop diagram this requires 

A'<3xKr n , (30) 

where we have taken m 3 / 2 ~ Mw- Clearly this bound is only approximate, since we do 
not know the precise values of the Yukawa couplings A and k, which we have taken here to 
be of 0(1). However, it should also be borne in mind that one would like to have control 
over the scale of electroweak symmetry breaking. That is we do not wish the mass of the 
W to depend strongly on the (unknown) physics at the Planck scale, i.e. on A'. In order 
to achieve this, the above bound should be tightened even further. 

The bound in eq. (|30|) is clearly incompatible with that in eq. (|(|) required for successful 
nucleosynthesis, and we conclude that the NMSSM at least in the models with 'minimal' 
Kahler potentials has either a domain wall problem or a hierarchy problem. 



6 A Solution to the Hierarchy Problem 

Is it possible that we can solve these problems by allowing the hidden and visible sectors 
to mix? In this section we shall see that the answer is yes for the destabilising divergences, 
but no for the naturalness problem. In other words, we are able to regain perturbative 
control over the scale of electroweak symmetry breaking, but we find, quite generally, that 
certain couplings must be set by hand initially to be small. This leads to a naturalness 
problem of at least one part in 10 9 . 

In order to eliminate destablising divergences, we must drop our insistence on mini- 
mality in the Kahler potential, by allowing the hidden and visible sectors to mix. In this 
case models similar to the NMSSM can be constructed. We use a mechanism similar to 
that in ref.0, and find that models with (standard model) singlets can have naturally 
large N 2 , N 3 and fi terms. 

The Giudice-Masiero mechanism || seeks to solve the fi problem for the MSSM by 
generating it via the Kahler potential. That is we have 

G = y'yj + zz^ + {^-z^H x H 2 + h.c.) + M^ln 
V Mp\ / 

where the yi fields belong to the visible sector, and the z singlet field belongs to the 
hidden sector. Q is Kahler invariant. The label 'hidden' is justified when we take the 
"flat" limit Mpi — > oo in the effective potential (keeping M s 2 /Mpi fixed), and find that 
the z field, which acquires a vev of O(Mpi), decouples from the visible sector, apart from 
inducing soft SUSY breaking terms and a /i term proportional to a, via gravitational 
couplings. These are all of (9(M s 2 /M P i), where M s is the aforementioned scale of SUSY 



f(z)+g(y) 



(31) 
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breaking in the hidden sector which we introduce by hand. Having introduced a new 
coupling between the visible and hidden sectors, we must invoke some symmetry which 
forbids other couplings as well as a coupling Mj>iHiH 2 in the superpotential. This could 
be a Peccei-Quinn symmetry, a discrete symmetry, or a gauged or global R symmetry. 
In addition the presence of a new symmetry rules out the simplest version of the Polonyi 
model (which in view of its severe cosmological problems |23| might not be such a bad 
thing). 

For the next-to-minimal choice of Kahler potential above, the terms in the scalar 
potential are 



^scalar = (Jiff + ml/flof + ^ [tf 9i + (A - 3)# {3) + (B - 2)to 3 / 2 /^1#2 + h.C. 

where are the trilinear terms of the superpotential, rescaled according to 

<? (3) = (exp(zz^/2M^)}g^\ 
Here g is the new low energy superpotential including the fx term 



9 = 9 



(3) 



fxH^, 



and m 3 /2 is the gravitino mass 



m 3/2 = (exp(zzy2M^)f^), 



(32) 



(33) 



(34) 



(35) 



where the vev of = M s 2 /Mpi is set by hand such that M s ~ 10 11 GeV. The fx term is 
given by 

/ i / -i- \ 

(36) 





/M P1 f z \ 


M = 


am \ f / 



Applying the constraint of vanishing cosmological constant, the authors of Ref. found 

B = (2A - 3)/(A - 3) 

\fi\ = \ma{A-3)/V3\, (37) 



where A is the universal trilinear scalar coupling, A = y3(z/Mpi). Now let us apply the 
same mechanism to a model with MSSM singlets, N. The most obvious extension is to 
choose the Kahler potential 



Q = y'yj + zz^ + 



M, 



" z^E x E 2 + ^z ] N 2 + h.c. ] + Mp, In 



pi 



Mi 



pi 



f(z)+g(y) 



(3? 



where, in this case, f(y) is the superpotential of the NMSSM. The hidden sector field, 
z, has the opposite charge to N under the Z 3 symmetry so that the full theory is Z 3 - 
invariant. In this case Z 3 is broken spontaneously at the Planck scale and the resulting 
domain walls are presumably removed during inflation. The low energy scalar potential 
is 



Scalar = ftS'+mW+m* tf 9i + (A - 3)g^ + (B - 2)mfxE x E 2 + (B - 2)mfi'N' 2 + h.c 



-3^(3) 



(39) 
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where 



g = + f iH 1 H 2 + fi'N 2 
ma(A — 3) 



ma'{A - 3) 



(40) 



Notice that the low energy model has generally far more terms in its low energy lagrangian 
than the NMSSM. The latter (and the Z 3 symmetry) is in fact recovered when we let 
a = a' = 0; thus we can break the Z 3 symmetry by as much or as little as we like. 

Although this model has removed the problem of destabilising divergences (it is now 
no longer possible to write down any of the divergent tadpole diagrams), it does not quite 
solve the naturalness problem (i.e. the presence of small couplings unprotected by any 
symmetry), since there is still the coupling zN which is allowed under the Z 3 symmetry, 
and which no other symmetry can forbid. These may be set to zero by hand and will stay 
zero by virtue of the nonrenormalization theorem. 

One might wonder if by somehow extending the Kahler potential it may be possible 
to exclude these terms. As we now show however, this is not the case, and no matter how 
complicated we make the Lagrangian, the naturalness problem associated with the absence 
of the zN couplings stays with us. Consider the most general supergravity Lagrangian, 
in which the only requirement we make is that the superpotential contains the terms 

Sg(y) = -^N a N b N c + A ttbc (0(#i# 2 UiV c , (41) 

where a,b,c are indices representing some symmetry group (discrete or otherwise), and 
the couplings are holomorphic function of the hidden sector fields, £ a = z a /Mp\. The 
breaking of Z 3 symmetry in the visible sector by operators of dimension-5, requires that 
we also include at least one of the operators, 



A ab (£,£)AyV ( 



b 



A?(£,0N a N b N^ 

A abcd (C)N a N b N c N d 
A abcd {0NaNb{HiH2)cd 

A ^ ){HiH2)ab{HiH2U (42) 

where the first four operators give dimension-5 operators if they appear in the Kahler 
potential or superpotential, but the last three operators must appear in the superpotential, 
hence their couplings are holomorphic functions of the hidden sector fields. If we make 
the assumption that the couplings k abc and X abc are invertible, then corresponding to each 
of the operators above, there is an additional invariant operator which is some function 
of the hidden sector fields multiplied by N a . These are, respectively, 

Ai b (k-y abc N c 
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Al^A- 1 )^ 
Kt{k- l ) abd {k- l )^ e N e 

A^k-xu^Uik-y^Ng 

A abcd {\- l ) abe {\~ l ) cdf {k- l )^Ng. (43) 

The least damage to the effective potential occurs if these terms appear in the Kahler 
potential, in which case we find terms of the form 

ml /2 Mpi<f) N + h.c. (44) 

appearing in the effective potential. Thus the natural scale of the singlet vev is ~ 10 11 
GeV and since it should be less than the electroweak scale, this constitutes a naturalness 
problem of at least one part in 10 9 . 



7 Conclusions 

Before we conclude there are a few escape clauses which should be mentioned none of 
which however are very appealing: 

1. The most obvious is to introduce the [i term into the superpotential by the mech- 
anism of Ref. || and simply set to zero all of the operators which might give N a 
large vev. (Although this appears to be rather unaesthetic, one might remark that 
the naturalness problem which results is no worse than that already with us due to 
the smallness of M s compared to the Planck mass. Since the "unnaturalness" is of 
the same order, it may even be possible to construct the Kahler potential so that 
the two naturalness problems are connected.) 

2. Alternatively one can invoke inflation at the weak scale to remove all the domain 
walls, just as has been suggested in the context of other unwanted relics, e.g. string 
moduli ||25|| . However such a scenario must be very finely tuned — the domain walls 



must be adequately diluted without erasing the density perturbations generated by 
inflation at the GUT scale p6fl ). (Although density perturbations are also generated 
during weak scale inflation, the small value of the Hubble parameter would make 
these too small to account for the microwave background anisotropies observed by 
COBE.) Secondly, the reheat temperature must be high enough for both successful 
baryogenesis and nucleosynthesis. We are not aware of any compelling candidate 
for the required scalar field. 

3. The Z 3 symmetry could be broken at a high scale, M contrived , in the visible sector and 
explicit Z% breaking terms induced. This is similar to the solution to the hierarchy 
problem discussed earlier, with the "advantage" that the fine tuning is driven by 
^contrived rather than Mp\. However this will still entail a fine tuning of approxi- 
mately one part in 10 12 , since in order for the walls to be inflated away, M contrived 
should exceed the inflationary scale of ~ 10 14 GeV as deduced from normalisation 



13 



to the COBE data [2B|. Otherwise one would have to invoke a second epoch of 
inflation at an intermediate scale, with its own attendant problems (see above). 



4. The Z% symmetry could be made anomalous by adding extra fields to the theory 
which couple to SU(3) C (for example an additional generation). In this case the 
symmetry is broken non-perturbatively at the quark-hadron phase transition, and 
the walls collapse very soon thereafter [pTf] . However, it is difficult to see how this 
constitutes a solution to fine-tuning, since at the same time it seems to preclude a 
solution to the strong CP problem as discussed in Ref . ||27[| . 

5. The Z3 symmetry could be embedded in a continuous gauge or global group which is 
broken at some high scale. This is the Lazarides-Shafi mechanism ||29|| , in which the 
apparent discrete symmetry is a subgroup of the centre of the continuous group. In 
this case only U(l), SU(3n) (where n is an integer) and E 6 are suitable candidates 
(see for example Ref.[Q). After the electroweak phase transition, one expects only 
a network of walls bounded by strings to form and then quickly collapse 



To summarize, we have shown that the domain wall problem in the NMSSM causes 
it to be ruled out on cosmological grounds unless we break the Z3 symmetry of the 
model explicitly. The breaking may be driven by terms which are non-renormalisable 
and have no direct effect on the low energy theory. However their introduction will in 
general generate terms which destabilise the hierarchy. In models with "minimal" Kahler 
potentials, we have shown that there are no non-renormalisable operators which can be 
added to the superpotential with a coefficient which is simultaneously large enough to 
solve the cosmological problem and small enough to avoid reintroducing the hierarchy 
problem. Furthermore, if any of these operators are allowed by the symmetries of the 
theory at the supergravity scale, then there is no possible symmetry which could prevent 
the existence of an operator zN in the superpotential whose coefficient must be £ 10~ 17 . 
If we allow mixing between the hidden and visible sectors, the reintroduction of the 
hierarchy problem can be avoided, and the naturalness problem can be formulated in a 
way very similar to the \x problem in the MSSM. However, even here we must arbitrarily 
select coefficients of certain dangerous operators to be of 0(1O~ 9 ) or less once we have 
aranged for a /i parameter of a reasonable size, and we have also reintroduced the /x 
problem which the model was at least partly designed to solve. Thus we conclude that 
the parameters in the NMSSM must be very strongly fine tuned if we wish to avoid both 
the cosmological problems associated with domain walls and the hierarchy problem, and 
hence that the model suffers from severe naturalness problems. 
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Figure Captions 
Figure la,b 

An example of a wall configuration with the singlet vev x = 10u. Here we have chosen tan (3 = 2, 
X = k = 0.2, A k = A x = 200GeV. Total surface energy density is 8.6 x 10 8 GeV 3 . Figure la 
shows the values of the three scalar fields as a function of position; Figure lb shows the energy 
density in the wall relative to the vacuum. 
Figure 2 

The function f(T) (see eq.12) related to friction plotted against temperature, as discussed in 
the text. 
Figure 3 

A typical example of the evolution of the wall network with pressure. The figure shows the 
wall network at four times separated by intervals of lCP 10 sec, with a pressure of term of order 
aMyy/Mpi, beginning at the time when pressure starts to dominate the evolution. 
Figure 4a,b,c 

The three dangerous diagrams for each of the three operators which can destabilise the hierarchy. 
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